Factors that predict better synchronizability on complex networks 
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While shorter characteristic path length has in general been believed to enhance synchronizability 
of a coupled oscillator system on a complex network, the suppressing tendency of the heterogeneity 
of the degree distribution, even for shorter characteristic path length, has also been reported. To 
see this, we investigate the effects of various factors such as the degree, characteristic path length, 
heterogeneity, and betweenness centrality on synchronization, and find a consistent trend between 
the synchronization and the betweenness centrality. The betweenness centrality is thus proposed as 
a good indicator for synchronizability. 
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In recent years, diverse systems in nature have been 
observed to exhibit the characteristics of complex net- 
works, drawing much attention to complex network sys- 
tems [l|, y| ■ Previous studies have mostly been focused 
on structural properties of the networks rather than dy- 
namical ones, even though dynamical properties are also 
very important for understandingthe systems on com- 
plex networks. In a few studies |J 1^, l^i on the other 
hand, the dynamical system of coupled oscillators has 
been considered on complex networks and collective syn- 
chronization displayed by the system has been investi- 
gated. There it has been found that shorter character- 
istic path length tends to enhance synchronization. In 
contrast to this, a recent paper 6], investigating the ef- 
fects of heterogeneity of the degree distribution on collec- 
tive synchronization, reported that synchronizability is 
suppressed as the degree distribution becomes more het- 
erogeneous, even for shorter characteristic path length. 
These different results then raise a question as to syn- 
chronization on complex networks: What is the most 
important ingredient for better synchronizability? 

As an attempt to give an answer to this, we in this 
paper consider a system of coupled limit-cycle oscillators 
on the Watts-Strogatz (WS) small- world network [3, and 
investigate collective synchronization of the system. We 
pay particular attention to how the synchronization is 
affected by various factors such as the maximum degree, 
characteristic path length, heterogeneity of the degree 
distribution, and betweenness centrality. Here the col- 
lective synchronization is explored via the eigenvalues of 
the coupling matrix, which describes the stability of the 
fully synchronized state 0, . 

The WS small-world network is constructed in the fol- 
lowing way j3|: We first consider a one-dimensional reg- 
ular network of N nodes under periodic boundary con- 
ditions, with only local connections of range r between 
the nodes. Next, each local link is visited once, and with 
the rewiring probability p it is removed and reconnected 
to a randomly chosen node. At each node of the small- 



world network built as above, an oscillator is placed; a 
link connecting two nodes represents coupling between 
the two oscillators at those two nodes. We now investi- 
gate the synchronization of the coupled oscillators on the 
small-world network with given r and p. Describing the 
state of the ith oscillator (i.e., the one at node i) by Xi, 
we begin with the set of equations of motion governing 
the dynamics of the N coupled oscillators: 

N 

X, = F{xi) +kY,M,jG{xj), (1) 

where Xi — F{xi) governs the dynamics of individual 
oscillators (i.e., with coupling strength K = {)) and G{xj) 
makes the output function. The N x N coupling matrix 
Mij is given by 

{ki for i = j 
-1 for J G A, (2) 
otherwise, 

which lacks the translational symmetry due to the pres- 
ence of shortcuts on the WS small-world network. In 
the case of a locally-coupled (hypercubic) network with 
the coordination number z (= 2D), the coupling matrix 
Mij has the value z on the diagonal and —1 on the z 
off-diagonals adjacent to the diagonal. 

The eigenvalues of the coupling matrix have been 
widely used to determine the linear stability of the fully 
synchronized state {xi = X2 = ■ ■ ■ = xn) Whereas 
the smallest eigenvalue, denoted by Aq, is always zero, 
the ratio of the maximum eigenvalue Amax to the small- 
est non-vanishing one Amin may be used as a measure for 
the stability of the synchronized state, with larger values 
of the ratio Apax/Amin corresponding to poor synchro- 
nizability |5|, t3 ■ For a general D-dimensional hypercubic 
network of linear size L, the eigenvalues are given by 

V.}=4f]sin2^, (3) 
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FIG. 1: Behavior of the ratio of the maximum eigenvalue 
Amax to the smallest non-vanishing eigenvalue Amin with the 
rewiring probability p. As p is raised, the ratio is shown to 
reduce, yielding better synchronization. 
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FIG. 2: Behavior of the variance crfc^ of the degree distri- 
bution with the rewiring probability p. Inset: characteristic 
path length £ vs. the rewiring probability p. Here and in 
subsequent figures lines are merely guides to eyes. 



where rua = 0, 1,-- - ,L—1 and = N. It is then 
obvious that the eigenvalue ratio behaves as Amax/ -^min ^ 
N^/° and grows large in the thermodynamic limit {N — > 
oo). It is thus concluded that the fully synchronized state 
is not stable in any D-dimensional regular network. 

The eigenvalue ratio for the WS network of nodes 
N = 2000 and range r = 3 is obtained numerically and 
its behavior with the rewiring probability is exhibited in 
Fig. n where the average has been taken over 100 differ- 
ent network realizations. The network size N has been 
varied from 100 to 2000, only to give no qualitative dif- 
ference. As the rewiring probability p is increased, the 
ratio A,nax/-^min IS obscrvcd to decrease, which implies 
enhancement of synchronizability. 

To explore how structural properties of the underlying 
network affect synchronization of the system, we now ex- 
amine such properties as the characteristic path length, 
the betweenness centrality, and the variance of the degree 
distribution, which is a measure of the heterogeneity of 
the degree distribution. Figure |21 displays the behavior 
of the variance = {N-^ h^) " ((^"^ T,^ hf) of 
the degree distribution for the WS small-world network 
with the same size N = 2000 and range r = 3 as that 
in Fig. ^ depending on the rewiring probability p. As p 
is increased, the variance Ok^ grows, which implies that 
the degree distribution becomes more broad and hetero- 
geneous and that nodes of larger degrees appear. In the 
inset of Fig. [3 the behavior of the characteristic path 
length [i| 
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(4) 



where (■ • • ) denotes the average over different realizations 
of the network and dij the length of the geodesic between 
i and j, is shown as a function of the rewiring probability 



p. As expected, the characteristic path length £ is ob- 
served to decrease as the heterogeneity of the degree dis- 
tribution (or the rewiring probability) is increased. The 
results shown in Figs. ^ and 12 imply that the synchro- 
nizability on the WS network is improved as the hetero- 
geneity of the degree distribution is increased or as the 
characteristic path length is decreased, which differs from 
the behavior observed in scale- free networks 0. 

The synchronizability has been shown to be related 
with the load or betweenness centrality on nodes [a. The 
betweenness centrality of node n is defined to be M 0, 

MM- 
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(5) 



where gij is the number of geodesic paths between nodes 
i and j and gi„j is the number of paths between i and 
j passing through node n. The summation is to be per- 
formed over all pairs of nodes such that i, j ^ n and 

To get an idea of the betweenness centrality, which 
measures how many geodesies pass through a given node, 
we first consider locally-coupled regular networks, for 
which the average betweenness centrality is given by [Tlj 
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(6) 



Among the N values of i?„'s, the maximum value B'^'^^ 
has been shown to be related with synchronizability 
although this close relation has not been stressed before. 
For a £)-dimensional local regular network, we have B = 
j^max ^YiQ characteristic path length £ ~ iV^/^, which 
yields 



B" 



(7) 
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FIG. 3: Behavior of the maximum value TJ™''" of the be- 
tweenness centrality with the rewiring probability p. 73™^" is 
shown to decrease as p is increased, corresponding to more 
heterogeneous degree distributions. 

i.e., B"^^'^ ~ N^, 7V3/2^ and TV^/s ^j. ^^le spatial dimen- 
sion D = 1, 2, and 3, respectively. Thus the maximum 
value B'^'^^ increases algebraically vsfith the size N of the 
regular network, although the exponent reduces with the 
spatial dimension D. 

Returning to the WS small-world network, we com- 
pute the betweenness centrality Bn via a modified ver- 
sion of the breadth- first search algorithm iHJi . We then 
obtain the maximum value B™'^^ at various values of the 
rewiring probability, and display the result for a network 
of nodes N = 2000 in Fig. |21 where the average has been 
taken over 100 different network realizations. The num- 
ber of nodes N has also been varied from N = 100 to 
N = 2000, which does not yield any qualitative differ- 
ence. Figure 121 shows that the maximum load on a node 
reduces as more shortcuts are introduced. In general, on 
a usual scale-free network, larger values of B"^'^^ corre- 
spond to larger values of the degree. For comparison, we 
also investigate the maximum degree k'^^^ on our small- 
world network and display in Fig. 21 its behavior with the 
rewiring probability p. The increase of k"^'^^ with p indi- 
cates the opposite trend between B"^^^ and k"^'^^, unlike 
the case of a scale- free network. Note also that Fig. Alto- 
gether with Fig. 121 implies the increase of the maximum 
degree with the heterogeneity of the degree distribution. 

The results shown in Figs.QtoEllead to the conclusion 
that synchronization on the WS network is enhanced as 
the heterogeneity of the degree distribution is increased, 
as the characteristic path length decreased, as the max- 
imum betweenness centrality decreased, or as the max- 
imum degree increased. Remarkably, the effects of the 
heterogeneity of the degree distribution as well as those 
of the characteristic path length (see Fig^^EJ differ from 
the results for other classes of networks "(f. Namely, in 
the case of the network studied here, larger heterogeneity 
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FIG. 4: The maximum degree fc™^^ versus the rewiring prob- 
ability p on a WS small-world network. It is shown that fc™™ 
increases as the degree distribution becomes more heteroge- 
neous. 



of the degree distribution or shorter characteristic path 
length does not improve synchronizability. On the other 
hand, the effects of the maximum betweenness central- 
ity B"^'^^ appear to be consistent with those in Ref. 
Synchronizability is always improved as B"^'^^ is reduced. 
Accordingly, the betweenness centrality is proposed as 
a suitable indicator for predicting synchronizability on 
complex networks. Regarding the maximum degree k"^^'^, 
Fig. 01 (together with Fig. ^ indicates that synchroniz- 
ability enhances with k"^^^. This behavior on the WS 
small-world network is also in contrast with that on a 
usual scale-free network, where B"^^^ and fc™^'' behave 
similarly , and accordingly, smaller values of /c™'*'^ are 
expected to give better synchronizability. 

Then why does the maximum value B''^^^ of the be- 
tweenness centrality strongly affect synchronizability on 
networks? An intuitive argument goes as follows: Sup- 
pose that two groups of highly linked nodes are con- 
nected via a few nodes. Among those a few connect- 
ing nodes, the information of the synchronized state 
passes through the node which has the maximum value 
of the betweenness centrality, B"^'^^. It tends to get over- 
loaded since most paths go through it, which in turn 
leads to loss of information. Accordingly, synchroniz- 
ability is expected to become reduced as the value B"^'^^ 
is increased. This argument is consistent with that of 
Ref. 0. To check this argument, we have examined 
synchronizability on the WS network before and after 
the removal of the node having B'^'^^, by means of the 
eigenvalue ratio Amax/Amin- In Fig- 5, the difference 

^ — (Aniax/Ajnin) after removal (Ainax/Ainin)bcforc removal is 

plotted against the rewiring probability p. Squares in 
Fig. 5 represent the difference S for the removal of the 
node having B'"^'^^, and indicate negative values regard- 
less of the rewiring probability p. Those negative val- 
ues imply that Amax/Amin is reduced after the removal 
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FIG. 5: Behavior of the difference S of the eigenvalue ratio 
is displayed with the rewiring probability p. The data (open 
circles and squares) have been obtained for the network with 
the number of nodes — 2000 and r = 3. Due to strong 
finite size effects at p = 0.1, we use the data (filled circle and 
square) for large system size of A'^ = 5000. See the text for 
the explanation of the data symbols. 



load on a node, is reduced, which is consistent with the 
recent result for various networks We have also nu- 
merically investigated the effects of the node of the max- 
imum betweenness centrality B™^^ on synchronizability, 
and found that the node of B™^'^ is highly related with 
the synchronizability of the system, which supports the 
main conclusion of this paper. It is thus concluded that 
among the important factors for better synchronization 
on complex networks is a small value of the maximum 
betweenness centrality, rather than short characteristic 
path length or large heterogeneity of the degree distribu- 
tion. 
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Nos. R14-2002-062-01000-0 (B.J.K.) and ROl-2002-000- 
00285-0 as well as by the Ministry of Education through 
the BK2I Project (M.Y.C.). Numerical works were 
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Ajou University. 



of the node, which in turn implies that synchronizabil- 
ity of the system is enhanced after the removal of the 
node. In sharp contrast, circles in Fig. 5, which repre- 
sent the difference 5 for the removal of a node other than 
the one with B'^'^^, indicate values near zero. This im- 
plies that the eigenvalue ratio remains almost the same 
when an arbitrary node is removed and accordingly that 
synchronizability of the system is not much affected by 
the removal of an arbitrary node. The result displayed in 
Fig. 5 thus manifests that the node with iJ™^,'^ is closely 
related with the synchronizability of the system whereas 
any other node is not substantially related. 

In conclusion, better synchronizability for the WS 
small-world network is induced as the heterogeneity of 
the degree distribution is increased, as the characteristic 
path length is decreased, as the maximum betweenness 
centrality is reduced, or as the maximum degree is raised. 
We have found that the effects of the characteristic path 
length and of the heterogeneity on synchronization in the 
WS small-world network are different from those in the 
networks considered in Rcf. 6J. These differences seem 
to be related with the presence of hub structures in the 
network, which is under further investigation. Our result 
implies that shorter characteristic path length or larger 
heterogeneity does not always enhance the synchroniz- 
ability of the coupled system on a network. On the other 
hand, it has been observed that synchronization is always 
enhanced as the betweenness centrality, measuring the 
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